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Abstract: We explore various aspects of the correspondence between dimer models 



and integrable systems recently introduced by Goncharov and Kenyon. Dimer models 
give rise to relativistic integrable systems that match those arising from 5d M — 1 gauge 
l~^ theories studied by Nekrasov. We apply the correspondence to dimer models associated 

to the Y p '° geometries, showing that they give rise to the relativistic generalization 
of the periodic Toda chain originally studied by Ruijsenaars. The correspondence 
reduces the calculation of all conserved charges to a straightforward combinatorial 
problem of enumerating non-intersecting paths in the dimer model. We show how 
the usual periodic Toda chain emerges in the non-relativistic limit and how the Lax 
operator corresponds to the Kasteleyn matrix of the dimer model. We discuss how 
the dimer models for general Y p,q manifolds give rise to other relativistic integrable 
systems, generalizing the periodic Toda chain and construct the integrable systems for 
general Y p,p explicitly. The impurities introduced in the construction of Y p,q quivers 
are identified with impurities in twisted sl(2) XXZ spin chains. Finally we discuss 
how the physical concept of higgsing a dimer model provides an efficient method for 
producing new integrable systems starting from known ones. We illustrate this idea by 
constructing the integrable systems for higgsings of Y 4,0 . 
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1. Introduction 

Integrable systems have a dense web of connections with gauge theories and string 
theory. Integrable systems appear in a variety of contexts including Seiberg-Witten 
theory |]J, recent connections to the vacua of supersymmetric theories [0, and the 
calculation of the spectrum of anomalous dimensions and scattering amplitudes Q 
in super- Yang-Mills, to name a few. 

Recently, Goncharov and Kenyon discovered an exciting correspondence between 
integrable systems and dimer models ||. According to their correspondence, every 
dimer model defines an integrable system, whose conserved charges can be systemati- 
cally calculated from perfect matchings. 



The correspondence sheds new light on integrable systems, with applications and 
implications yet to be investigated. It provides new perspectives on integrable systems 
that are naturally related to dimer models, such as 4d quiver gauge theories, D3-branes 
probing toric Calabi-Yau 3- folds ||, mirror symmetry [[/J and quantum Teichmiiller 
space §0. 

In this paper, we investigate various aspects of the correspondence from a physical 
perspective. This work is organized as follows. Section 2 briefly reviews the relation 
between dimer models and quiver theories, toric Calabi-Yaus and integrable systems. 
Section 3 discusses the connection between integrable systems and 5d H = 1 and 4d 
N = 2 gauge theories. In Section 4, we apply the correspondence, constructing a 
relativistic generalization of the periodic Toda chain from the dimer models associated 
to Y p '° manifolds. We study the relation between the Kasteleyn and Lax operators, the 
non-relativistic limit of the integrable system and how dimer models for general Y p,q 
geometries produce alternative relativistic generalizations of the periodic Toda chain. 
We also discuss the connection to twisted si(2) XXZ spin chains with impurities. In 
Section 5, we introduce a practical method for generating new integrable systems based 
in higgsing and illustrate the method with explicit examples. We conclude and mention 
future directions in Section 6. 



2. Some Background 

In this section we provide a lightning review of various concepts used throughout this 
paper. When necessary, we indicate references for more thorough explanations. 

Dimers and Quivers 

Brane tilings, to which we will also refer to as dimer models, are bipartite graphs 
embedded in a two-torus. The dual of a brane tiling is a planar, periodic quiver. There 
is a one-to-one correspondence between brane tilings and periodic quivers |6j that is 
summarized in the following dictionary: 



Gauge Theory 




Brane Tiling 


gauge group 


-H- 


face 


chiral superfield 


-H- 


edge 


superpotential term 


-B- 


node 



Every term in the superpotential of the gauge theory is encoded in an oriented 
plaquette of the periodic quiver. Figure |l| exemplifies the correspondence for phase I 
ofFoH. 




node = W term 



face = gauge group 



edge = bifundamental 




Figure 1: Correspondence between the brane tiling and the periodic quiver for phase I of 
F . 



Dimers and Calabi-Yau 3-folds 

Quiver gauge theories that are described by brane tilings arise on the worldvolume of 
stacks of D3-branes probing singular, toric Calabi-Yau (CY) 3-folds. The CY geome- 
try emerges as the moduli space of vacua of the quiver gauge theory. The connection 
between dimer models and quivers has trivialized the determination of the correspond- 
ing CY geometry. GLSM fields in the toric description of the CY are in one-to-one 
correspondence with perfect matchings of the dimer model. As a result, points in the 
toric diagram correspond to (sets of) perfect matchings. This correspondence reduces 
the task of finding the CY geometry to computing the determinant of the Kasteleyn 
matrix H. 



Dimers and Integrable Systems 

The dynamical variables of the integrable system correpond to oriented loops in the 
brane tiling. One basis for such loops is given by the cycles going clockwise around 
each face w t (i = 1, . . . , N g , with N g the number of gauge groups in the quiver) and 
the cycles Z\ and Z2 wrapping the two directions of the 2-torus. 1,2 



1, one of the Wj's is redundant. This identity can also be exploited for simplifying 



1 Since njl'i t 
expressions. 

2 The analysis of some models, such as the ones in Section |[ can be considerably simplified by 
choosing a different basis. 



The Poisson brackets between basis cycles are 

{wi,Wj} = e Wi>Wj WiWj 

{z 1 ,z 2 } = l + e ZuZ2 (2.1) 

where e Xjy is the number of edges on which the loops x and y overlap, counted with 
orientation. Then, e WijW , is simply the antisymmetric oriented incidence matrix of the 
quiver. 

The classical integrable system can be quantized replacing the Poisson brackets by 
a q-deformed algebra of the form 

XiX^qte^XjXi, (2.2) 

where Xi = e Xl and q = e~ l27rh . 

Every perfect matching corresponds to a point in the toric diagram and defines 
a closed loop by subtraction of a reference perfect matching. This loop can then be 
expressed in terms of the basic cycles. When multiple perfect matchings correspond to 
a given toric diagram point, their contributions must be added. Goncharov and Kenyon 
showed that the commutators defined by ( |2.2|) and ( p.l|) result in a (0 + l)-dimensional 
quantum integrable system in which the conserved charges are given by: 



Casimirs: they commute with everything. They are defined as the ratio between 
contributions of consecutive points on the boundary of the toric diagram. 

Hamiltonians: they commute with each other and correspond to internal points 
in the toric diagram. 



In this paper we will not discuss the choice of reference perfect matching in detail, 
we instead refer the interested reader to . Different choices of the reference perfect 
matching correspond to shifts in the toric diagram. These overall shifts do not affect 
Casimirs (since they are defined as ratios of points in the toric diagram) but they 
modify the Hamiltonian(s). Models with zero or one internal point are insensitive to 
this choice, since they have zero or one Hamiltonians. The choice of reference perfect 
matching becomes important for models with more than one internal point and can be 
straightforwardly determined by demanding the Hamiltonians to commute. 

Following |II| (see also [|ll|] for applications), we define the magnetic flux through 



a loop, 7, in terms of edges in the tiling as 



fc-1 



X(w i; bi) 



"(7) = 11 w w .,, hV ( 2 -3) 



^ X(w m , ^ 



where the product runs over the contour 7 and b« and Wj denote black and white 
nodes. With this definition, the elements in our basis are Wj = v(j Wj ), Z\ = v(-y Zl ), 
and z 2 = f (7 Z2 ). Another brief but more complete summary of the work in || can be 
found in M. 



3. Integrable Systems from Dimers, 5d and 4d 

5d Gauge Theories and Dimers 

The integrable systems we are discussing can be derived from either dimer models or 
5d gauge theories. The main object underlying all constructions is the spectral curve 
£. 

Dimer models encode the quiver gauge theory (and also the geometry) on D3-branes 
probing a singular, toric Calabi-Yau 3-fold X. The toric singularity has a characteristic 
polynomial P(zi,z 2 ) = J2 a m,n2 z™ 1 z 2 2 , where (711,712) runs over points in the toric 
diagram. The mirror Calabi-Yau is given by P(z\, z 2 ) = W, W = uv. In this case, E 
corresponds to the Riemann surface sitting at W = in the mirror. 

Nekrasov |T2[ proposed a non-perturbative solution for 5d gauge theories compact- 



ified on a circle using relativistic integrable systems. From the Seiberg-Witten solution 
fll3"l, |i"4fl a connection to integrable systems was proposed in |j]. In particular, the 
spectral curve of the integrable system matches the Seiberg-Witten curve. Nekrasov's 
insight was to generalize the integrable system to a relativistic integrable system in 
order to determine the corresponding Seiberg-Witten curve for the 5d gauge theory. 
There are two possible ways to engineer such compactified 5d theories and £ plays a 
prominent role in both of them. First, we can construct them by wrapping an M5- 
brane on E. This M5-brane is a de-singularization of a web of (p, q) 5-branes in Type 
IIB, obtained after compactifying in a circle of radius /3 to pass to Type IIA and then 
lifting to M-theory |T5[ . Alternatively, these theories can be obtained in the low energy 



limit of M-theory on a Calabi-Yau 3-fold with vanishing cycles to decouple gravity 
|I6|, IT?] , [T~8| , jl9| . The Calabi-Yau is the same X of the dimer construction, so we see 
that E is also relevant from this perspective. Particle-like states arise from M2-branes 
wrapped around SUSY 2-cycles. 



4d Gauge Theories and the Non-Relativistic Limit 

The integrable systems constructed using dimer models (equivalently 5d gauge theories 
compactified on a circle) are naturally relativistic. This is reflected, for example, in the 
exponential dependence on momenta of the conserved charges. This fact is manifest in 
the form of Poisson brackets fl2.ip . 

The radius, (3, of the circle on which the 5d gauge theory is compactified plays 
the role of the inverse speed of light. The non-relativistic limit corresponds to taking 
(3 — > 0. 3 In this limit, we obtain a 4d, M = 2 gauge theory whose Seiberg-Witten 
curve is the spectral curve of the non-relativistic integrable system, which we denote 
a. 4 Figure |2| summarizes the connection between 5d gauge theories, dimer models, 
integrable systems and their non-relativistic (4d) limits. 



5d N=l gauge theory 
onS 1 

M5-brane wrapped on E 
M-theory on CY3 



4d N=2 gauge theory 

Seiberg-Witten curve a 



Relativistic 
Integrable System 

Spectral curve £ 



I 



pi— 



Non-Relativistic 
Integrable System 

Spectral curve a 



Dimer Model 

L inside the mirror 



Figure 2: Dimer models and 5d, TV = 1 gauge theories compactified on an S 1 correspond 
to different perspectives on the spectral curve X. The 4d, J\f = 2 limit of the gauge theory 
corresponds to a non-relativistic limit of the integrable system. 



4. The Periodic Toda Chain 

In this section we show the correspondence from || at work in explicit examples, using 



3 Also interesting is the decompactification limit f3 — » oo, in which the 5d perturbative solution [[20| 
is recovered 

4 As explained in Section 0, dimer models are in one-to-one correspondence with 4d (generically 
M = 1) quiver theories. These quivers thus provide other 4d gauge theories naturally associated to 
the integrable systems. In what follows, whether we refer to the 4d TV = 2 or quiver gauge theories 
should be clear from the context. 



dimer models to construct (relativistic generalizations) of the periodic Toda chain. The 
relativistic, periodic Toda chain was first introduced in plj and studied in connection 



with 5d gauge theories in ||12| . These integrable systems arise from dimer models, 
equivalently quiver gauge theories, associated to Y p ' q manifolds [p^j . The associated 
spectral curves, E, correspond to 5d, N = 1, SU(p) gauge theory with no flavors 
and different values of a quantized parameter controlling the cubic couplings in the 
prepotential. 

4.1 Y p '° Integrable Systems 

Integrable system 

Let us now consider general Y p, ° geometries. Our goal is to make contact with the 
relativistic Toda chain for arbitrary p. For concreteness, let us focus on the case in 
which p is even. The shape of the unit cell depends on whether p is even (rectangle) 
or odd (rhombus). The dimer model for the conifold corresponds to a square lattice 
p3| . The cone over Y p '° is a Z p orbifold of the conifold. As a result, its dimer model 
is given by a square lattice with an enlarged unit cell, given by p copies of the one for 
the conifold, as shown in Figure |^. 

The reference perfect matching and z\ and Z2 paths for Y p, ° are shown in Figure |j. 
Figure || shows the toric diagram for Y p, °, with even p. The system has a Z2 symmetry 
corresponding to choosing the opposite corner of the toric diagram as the reference 
perfect matching. The Z 2 symmetry interchanges the Hamiltonian with its dual |26| . 



The construction of the associated integrable system is considerably simplified by 
an appropriate choice of basis of 2p + 2 cycles, instead of Wi and Zj. Figure | shows 2p 
of the loops. There are two additional cycles, with winding numbers (—p/2 — 1, 1) and 
(— p/2 — 1, —1) around the z\ and Z2 directions, which correspond to the green points 
in Figure |5|. They are mapped to two Casimirs and are hence fixed. 

The non-vanishing Poisson brackets are 

{c k ,d k } = c k d k {c k ,d k+1 } = -c k d k {c k ,c k+ i} = -c k c k+1 (4.1) 

The Hamiltonian corresponds to the (—1,0) point in the toric diagram and is given by 

p 



H 1 = Y, ( ,a + d i ). (4.2) 



i=l 



Equations (|4.1|) and fl4.2|) precisely agree with the Poisson brackets and Hamiltonian 



for the general periodic relativistic Toda chain |25| p6| . The q and di variables can be 
expressed in terms of position and momentum variables with canonical commutation 
relations as follows 




U-U 




Hi Hi 

Figure 3: Brane tiling for Y p '° with even p. We indicate the numbering of nodes in blue. 



Ci = exp(pi -qi + q i+ i) 

di = exppi. (4.3) 

Using the q and di cycles, determining the additional (p — 2) higher Hamiltonians 
reduces to a straightforward combinatorial problem. The H n Hamiltonian, associated 
with the (—n, 0) point in the toric diagram, corresponds to the sum of all possible 
combinations of n of these cycles with the condition that they do not overlap or touch 
at any vertex of the tiling. For example, for Y 4 ' we have: 

H 2 = (cic 3 + c 2 c^) + (cx(ds + d^ + c 2 (d 4: + d\) + c 3 {d\ + d 2 ) + c A (d 2 + d 3 )) (4.4) 
+ (did 2 + did 3 + d\d± + d 2 d 3 + d^d^) , 



H 3 = (d 2 d 3 d^ + did 3 di + d\d 2 d^ + d\d 2 d 3 ) + (cid^ + c 2 d±d\ + c 3 d\d 2 + c^d 2 d 3 ) . 



■ti 



HP-*- 



ref. p.m. z\ Ti 

Figure 4: Reference perfect matching and z\ and Z2 paths for Y p '°. 




p/2+1 



Figure 5: Toric diagram for Y p, ° (shown in the figure for p = 6). The reference perfect 
matching is circled in red. By construction, its position in the (2:1,2:2) plane is (0,0). The 
green dots correspond to cycles with windings (— p/2 — 1, 1) and (— p/2 — 1, —1). 



4.2 Kasteleyn matrix 

Let us construct the Kasteleyn matrix for the dimer model in Figure [| It is convenient 
to label edges according to whether they are horizontal (H and H) or vertical (V and 
V). H edges are those at the center of the tiling and H are those crossing the edge of 
the unit cell. In addition: 



V : vertical with black node at top endpoint 

V : vertical with white node at top endpoint 

Finally, let us call Hi and Hi the edges connecting nodes i and n + i, Vi = Vi+i jP +i 
and Vi = Vi tP+ i + i. Subindices indicate nodes on the tiling and are identified mod(2p). 
In these variables, the Kasteleyn matrix is 



d, 
i=l,..,P 



Ci-i 
even i 



c, 
even i 



Figure 6: 2p of the cycles in a convenient basis for Y p, °. Notice that the c cycles only exist 
for even i. 



K n 
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p+1 p+2 


p + 3 ■ 


■■ 2p - 1 




2p 
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-H 1 -H 1 z 1 V x 

V 1 H 2 + H 2 z^ 1 


V-2 






V p z 2 


3 


v 2 










p — 1 
\ P 


Vp z 2 




Vr-x 


H p 


Vp-i 
+ H pZl -\ 



(4.5) 
Notice the alternating overall sign of the terms on the diagonal. 



4.2.1 Non-Relativistic Limit 



The Kasteleyn matrix (|4.5|) and the Lax operator of the non-relativistic periodic Toda 
chain |27|, ^] are strikingly similar. We now make this connection explicit by first 
expressing the edge variables in terms of coordinates and momenta and then taking the 
non-relativistic (i.e. small momentum) limit. 5 We define 



5 The definitions in (4.6) are reasonable, rather symmetric and simple. We later check that they 
arc indeed consistent with the w, commutation relations. 



10 



Vi 
H 



V,, = e 



qt-qi+i 



Hi 



= e (-iy P i/2 



(4.6) 



and 



z ± = e z z 2 = w (4.7) 

In the next section, we explain how these definitions and canonical Poisson brackets 
{Pii Qj} = $ij are i n agreement with the Goncharov-Kenyon Poisson brackets. 6 

Taking the small momenta limit (i.e. linear order in pi and z) of the Kasteleyn 
matrix ( |4.5|) , we conclude that 



with 



L p (w) 



det K p = det (L p (w) — z) 



pi e qi-Q2 

o<?i-<?2 p e Q2-qz 



oQ2-Q3 



(U 



e q P -qi w \ 



V 



(4.9) 



e Qp~i-<lp 

e q P -qi w -i e q P ~i-q P p p I 

which is precisely the Lax operator of the non-relativistic periodic Toda chain. Notice 
that the previous analysis nicely associates coordinates and momenta with the vertical 
and horizontal directions of the square lattice, respectively. 

4.2.2 Loop Poisson Brackets from Edges 

In the previous section we have expressed edges of the brane tiling in term of coordinates 
and momenta in such a way that the Lax operator of the periodic Toda chain is obtained 
from the Kasteleyn matrix of the dimer model by taking the non-relativistic limit. We 
now show how the commutation relations among loop variables given by the Goncharov- 
Kenyon rules are recovered from our edge definitions and the {pi, qj} Poisson brackets. 
Consider Y p '° with even p. Using equation (|2.3|) , we have 



odd i: 

even i: 



Wi 



H t Vr l H l+l Vr l 
H^ Vl H7 + \V t 



w 



w 



p+i 
p+i 



Hi Vr^H i+1 Vr l 



(4.10) 



l,...,p. From ( [L10D and the canonical Poisson brackets {pi,qj} = Sij (which 



become \pi,qj] = —ih5ij in the quantum theory) we can calculate 



^Changing from the basis of coordinates and momenta in (4.6) to the one in (|4.3|) is straightforward. 
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Wi w 2 w p 
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W 2 
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W p -! 
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to to 

to to 

to 
to to 
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-4 

-4 
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Wp+\ 

Wp+2 
Wp+3 

W 2 p-1 
V W 2 p 
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-4 
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4 

-4 


to to 

to to 
to 

to to 
to 

to to 



(4.11) 



which, modulo an unimportant overall scaling, are exactly the Poisson brackets that 
follow from Goncharov-Kenyon prescription! This scaling can be absorbed in the value 
of h in the quantum theory. It is important to notice that the final result depends 
crucially on the details of ( |4.6| ), such as the sign of the exponents, which are also vital 
for obtaining the correct non-relativistic limit. An interesting example that depends 
on these details is {iVi,w p +i+i} = 0. 7 



4.3 More Relativistic Generalizations of Toda from Y p,q 

The integrable models based on Y p, ° are not the only possible relativistic generalizations 
of the periodic Toda chain. In fact, as we now discuss, all Y p,q geometries give rise to 
valid generalizations. The corresponding toric diagram is shown Figure |7|. 8 

M-theory on the CY cones over Y p,q geometries gives rise to 5d, J\f — 1, pure SU(p) 
gauge theories. The theories differ in the value of c c i, which is a quantized parameter of 
the theory that controls the cubic couplings in the exact quantum prepotential, related 
to a five dimensional Chern-Simons term 



19 . For Y p,q , we have c r 



Q- 



The distinction between theories with different values of c c i disappears when taking 
the non-relativistic limit. This fact is clearly manifest at the level of the spectral curve. 



7 Of course, the variables in (L6) also reproduce the Poisson brackets ( [4.l[ ) for the basis of cycles 



we considered in Section \i.l\ Those cycles are typically written in terms of coordinates and momenta 
such that horizontal lines are equal to e Pi and squares are equal to e 9 * -9 ^ 1 . The pi and qi considered 



in this section are the ones in which the Lax operator takes the form (4.9), but can be mapped to p\ 
and q[ by an appropriate change of variables. 

8 The toric diagram in Figure @ can be put into this form by an SL(2, Z) transformation. 
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(-l.p-q) 




(0,0) (1,0) 

Figure 7: Toric diagram for the cone over Y p ' q . 



The spectral curve associated with Y p,q and the manipulations for taking its non- 
relativistic limit have repeatedly appeared in the literature (see for example |3D| , which 
we now follow). 9 The spectral curve, E, can be written as 

p—q 

z x + a^— + P p (z 2 ) = (4.12) 

with P p {z2) a degree p polynomial. We can rewrite E as 

v 
y 2 = Y[{z 2 - e^f - Ae- tB z p 2 ~ q (4.13) 

i=i 
In order to take the non-relativistic (4d) limit, we define 

z 2 = e Px , ef = e^^ 1 , e~ tB = ( — j (4.14) 



2 



■ 



with a^j+i = a.i — Oj+i. From our perspective, it is clear that not only ^2 but also the 4>i 
variables and ts must be rescaled when taking the non-relativistic limit since they are 
controlled by the Wi variables, which are momentum dependent. In the — > limit, E 
becomes 

.2 



V 



p f\\ 2p 

n^-«M+i) 2 -4 (±\ , (4.15) 

i=\ ^ ' 

which is the Seiberg-Witten curve for the pure M = 2 SU(N) gauge theory. As we 
have anticipated, the dependence on q has disappeared. Reversing the reasoning, we 
conclude that all Y p,q manifolds (i.e. for arbitrary values of q) give rise to integrable 
systems that can be considered relativistic generalizations of the periodic Toda chain. In 
the next section we discuss the Y p,p geometry and its corresponding integrable system. 



9 In fact, this correspondence of the spectral curve is the original reason why we identify Y p,q dimcr 
models as giving rise to relativistic generalizations of the periodic Toda chain. 
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4.4 Y p ' p Integrable Systems 



We now construct a new infinite family of relativistic integrable systems associated 
with Y p,p . From the discussion in section [4.3| , we know that these systems also reduce 
to the periodic Toda chain in the non-relativistic limit. 

The cone over Y p,p is the C 3 /Z 2p orbifold. Its brane tiling consist of two columns 
of p hexagons as shown in Figure & 





Figure 8: Brane tiling for Y p ' p . A convenient basis for Y P)P is given by the Ui (up) and d% 
(down) cycles, i = 1, . . . ,p. 



Figure [5] shows the toric diagram for Y p,p , where the reference perfect matching is 
indicated with a red circle. As for Y p '°, the analysis of these models is simplified by 
a convenient choice of basis for closed cycles. Figure [S] shows 2p of them. There are 
two additional cycles with windings (— p, 1) and (— p, —1) along the (z\, z-i) directions. 
They correspond to the green points in Figure ^§ and are fixed by the Casimirs. 

The Hamiltonian associated to the (—1,0) point in the toric diagram takes the 
simple form 



Hi = ^{ui + di 



(4.16) 



i=\ 



Similar to our analysis of Y p '°, determining the Hamiltonians for the (— n, 0) points is 
very simple in the U\ and di basis. Finding the n th Hamiltonian reduces to determining 
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Figure 9: Toric diagram for Y p ' p (shown in the figure, p = 3). The reference perfect matching 
is circled in red. By construction, its position in the (z±,Z2) plane is (0,0). The green dots 
correspond to cycles with windings (— p, 1) and (— p, —1). 



all possible combinations of n cycles that do not overlap or intersect at nodes of the 
tiling. For example, the higher Hamiltonians for Y > are 



H 2 = U\U 2 + U1M3 + M1U4 + U 2 U 3 + U 2 U A + ^3^4 

+ d\d 2 + d\d 3 + didi + d 2 d 3 + d 2 d A + d^d^ 

+ U\d-3 + itidi + u 2 d4 + u 2 d\ + u^di + u 3 d 2 + u^d 2 + u^d 3 , 

#3 = M1U2U3 + U\U 2 U A + U1M3M4 + «2^3M4 

+ U\d^d^ + u 2 didi + u 3 d\d 2 + Uid 2 d 3 
+ d\d 2 d 3 + d\d 2 d^ + did^d^ + d^^d^ 
+ d x u 2 U3 + d 2 u 3 u 4 + d^u^ui + d A uiU 2 , 

i?4 = M1U2W3M4 + d\d 2 d 3 d^ . (4-17) 

4.5 F p ' g Integrable Systems as Spin Chains 

We have identified the family of integrable systems associated to the F 1 ''' 7 dimer models 
in the previous section. Previously the integrable systems identified with the Y p,q 



spectral curves were described as twisted sl(2) XXZ spin chains with impurities fl3T |. 
We will now explain the equivalence of these apparently different descriptions. 

The sl(2) spin chains are described by N "spins" \l/j where i — 1, . . . , N. The spin 
operators satisfy the commutation relations 

{S±, So} = ±5±, {S+, S-} = sinh25'o 

where the raising and lowering operators S± — Si ± iS 2 are denned as usual. Integra- 
bility of the spin chain can be shown starting from the auxiliary linear problem for the 
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Lax matrix 

with twisted boundary conditions implemented by the identification 

The spectral curve of the spin chain is given by the determinant of the transfer matrix, 

det(T(A) + wl) = 0, 

where the transfer matrix is defined by the product of the two-by-two Lax matrices 

T{X) = L N {X)...L 1 (X). 

Impurities are added to the spin chain by performing a site-dependent shift of 
variables for the Lax matrices Lj(fi). Y p ' q quiver gauge theories can be constructed 



by starting from Y p ' p and adding (p — q) impurities [ 22| . Amusingly, the spin chain 
and quiver impurities are precisely the same. This reflects the well-known phenomena 
that the same integrable system can have different Lax representations. Since the Lax 
matrices for all of the Y p,q quiver gauge theories are tridiagonal 0, we can re- write the 
determinant of the Lax matrix in spin-chain form using the following identity 



det 



(ax bi c \ 

ci '•• '•• 

' • • ' • • b N -i 
\b c N _ x a N J 



{-l) N ~ l ]J b 3 + Y[ cA + TrL N L N ^ ...L x 



where 

'a, -bj-iCj- 

v l C 

Thus we can re-write the the spectral curve as 



■' - [ 1 



det(T(A) + wl) = w 2 + wTiT(X) + det T(A). 
Under a suitable change of variables, it should be straightforward to show that this 



representation matches the XXZ form proposed in |31|| . For the relativistic, periodic 
Toda chain this change of variables appears in |29| . 
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5. Generating New Integrable Systems via Partial Resolution 

In this section we explain how to determine the integrable system associated with a 
partial resolution, given the one for the parent theory. This is a useful way of obtaining 
new integrable systems from known ones. Starting from a relatively complicated ex- 
ample, partial resolution provides a practical way of deriving new integrable systems, 
much faster in practice that going through the process of expressing the new system in 
terms of loop variables. This method is so efficient that it is natural to expect that it 
has a counterpart in the integrable system literature. 

We will focus on minimal partial resolutions, which correspond to removing ex- 
tremal perfect matchings (i.e those located at corners of the toric diagram) one at 
a time. 10 In addition, it might be necessary to simultaneously remove non-extremal 
perfect matchings. The result of this process is a new toric diagram in which the 
multiplicity of each new extremal perfect matching is one and some of the internal 
multiplicities might also change. From a quiver point of view, partial resolution cor- 
responds to turning on non-zero vacuum expectation values (vevs) for bifundamental 
fields Xij and then higgsing nodes % and j. 11 This operation might result in mass terms 
for some fields, which can be integrated out using their equations of motion. From 
a dimer model perspective, partial resolution corresponds to removing from the tiling 
the edges associated to the fields with non-zero vevs. As a result, some of the adjacent 
faces in the tiling are merged into a single new face. The integration of massive fields 
maps to the removal of 2-valent vertices that might be generated in the process by 
condensing the nodes at the endpoints of the two edges terminating in them. 

The standard understanding of partial resolutions using dimer models is that all 
perfect matchings containing the edge associated with X^ are removed M. Following 
the connection between bifundamental fields and perfect matchings given in , we see 
that all these perfect matchings (which are interpreted as GLSM fields) need to acquire 
a non-zero vev in order for the bifundamental to get one. 

In the integrable systems context, the fundamental objects are closed loops rather 
than individual perfect matchings. In fact, even the reference perfect matching might 
disappear when applying the discussion in the previous paragraph. It is straightforward 
to adapt the previous reasoning to loops: partial resolution removes all loops that 
contain an edge that gets a non-zero vev. 

From a brane tiling perspective, the edge associated to X^ is deleted and faces i 
and j are combined into a single face. Consequently, we start from two cycles Wi and 



10 The implementation of more general partial resolutions using dimer models has been discussed in 
great detail in |32]. We will restrict ourselves to minimal ones in this paper. 

n Notice that, given a quiver, not all possible higgsings correspond to consistent partial resolutions. 
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Wj and end with a combined cycle Wi/j = WiWj as shown in Figure [K 




o^o 




*v* 




•V* 



Figure 10: Combination of gauge group cycles when higgsing by a vev associated to the 
blue edge. 



The following rules produce the integrable system for the partially resolved geom- 



etry: 



1) Remove loops that contain an edge with a non-zero vev. 



2) Re-express the surviving loops with the replacement (WiWj) — > Wi/j. 



More practically, the loops that are removed in rule (1) are those that cannot be 
re- written using rule (2). 

In some cases, a Z{ path can involve an edge that is removed when higgsing, as 



shown in Figure |TT| If so, the path can be redefined by using a Wj to make the path 
wiggle appropriately, avoiding this edge. Equivalently, we could have chosen a different 
set of Zi paths in the parent theory such that they do not involve higgsed edges. 



o- 



3^0 




V* 




Figure 11: Redefinition of one of the paths that winds around the T 2 , z% = ZiWj, after partial 
resolution by turning on a vev for the blue edge. 



5.1 Examples: Partial Resolutions of Y 



•1,0 



We now illustrate our ideas in the explicit case of partial resolutions of Y 4 ' . Partial 
resolution can either preserve or reduce the genus of the spectral curve (i.e. the number 
of Hamiltonians) . The examples that follow exhibit the latter behavior. 

In this section we depart from the notation for edges we used for general Y p '° in 
Section |4.2| , which was specially devised for giving the Kasteleyn matrix a nice form. 
Here our emphasis is on higgsing, so we explicitly indicate the gauge groups under 
which bifundamentals are charged using subindices. Vy and V^ indicate vertical edges 
in the first and second columns of the brane tiling, respectively. Horizontal edges are 
denoted Hij. 



■■4- 



-0- 



-6- 



V, 



8 



-0- 



V — : Hi 



Figure 12: Brane tiling for Y 



1.0 



Figure [13] shows the resolutions we will consider. The number of gauge groups in 
the associated quivers is given by twice the area of the toric diagrams. This implies 
that the number of bifundamental expectation values that need to be turned on is equal 
to twice the decrease in area of the toric diagram. 

The starting point is the integrable system for Y A '°. In the table below, (711,712) 



gives the 



' ni 712 > 

/C-\ /in , 



contribution. Every term in a given contribution arises from a loop in 



the tiling or equivalently from a perfect matching. Hamiltonians correspond to internal 
points and Casimirs are given by the ratio of consecutive points on the boundary. 



Instead of the basis of cycles used in Section |Q| , here we use the W{ (which correspond 
to gauge groups), z\ and z 2 basis, since it makes higgsing more transparent. 
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V15, V37, V62, V84 

Figure 13: Toric diagrams for Y 4,0 and various of its partial resolutions. We indicate the 
multiplicity associated with each point and the non-zero vevs that need to be turned on. 
Higgsing takes the number of Hamiltonians from 3 in the original theory to a) 2, b) 1 and c) 
0. 



(ni,7i2) 


Loops 


(0,0) 


1 


(-1,0) 


W4 + W4W8 + W4W7WS + W3W4W7W8 + W2W3W4W7W8 
+W2W3W4WQW7W8 + W2W3W4W5WQW7W8 + 1 


(-2,0) 


Wi l W§ W4 + W4W8 + W^W^Ws + W^ 1 W^ 1 W4W 8 

+Wi 1 W^ 1 Wq 1 W4Ws + W4W7W8 + W^WiWjWs + Wi W5 W4W7W8 

+W3W4W7W8 + Wi 1 W 3 W4W 7 W8 + W^ 1 W^ 1 W 3 W4W 7 W 8 + W2W3W4W7W8 

+w± 1 w^ 1 + w 3 wlw 7 w 8 + w 2 w 3 wlw 7 w s + w 3 w\w 7 wl 


(-3,0) 


Wi 1 W$ 1 W4W8 + Wi 1 W^ 1 Wq 1 W4W8 + Wi W^ 1 W4W 7 Ws + Wi 1 W^ 1 W 3 WaW 7 W8 
+Wi 1 W^ 1 W 3 wlw 7 W8 + W 3 w\w 7 w1 + Wi WaW^WflVg + Wi W5 W3W4W7WI 


(-4,0) 


Wi 1 W^ 1 W3wlw 7 wl 


(-2,1) 


W^ l W4W 7 W8 


(-2,-1) 


W 2 W 3 wlw 7 W8 



(5.1) 
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Higgsing a 

Model (a) in Figure is obtained by giving non-zero vevs to both H 23 and H^q. There 



are other choices of expectation values that lead to the same result. Figure |14] shows 
the corresponding brane tiling. 



■■4- 



-0- 



2/3 Q 



0- 



5/6 



-0- 



Figure 14: Brane tiling for model (a) obtained by higgsing Y 



1.0 



The resulting integrable model is given by: 



(ni,n 2 ) 


Loops 


(0,0) 


1 


(-1,0) 


ti?4 + W4W8 + W4W7W8 + W2/3W4W7W8 
+W 2 /3W 4 W 5 /gW 7 W 8 + 1 


(-2,0) 


W4W8 + w^WiWs + w^w^hwiWg + W4W7W8 

+Wi 1 W i W 7 W 8 + W 2 /3W4W 7 Wg + W 2 /3w\w 7 Ws 


(-3,0) 


W^W^W^Ws 


(-2,1) 


W^W^WiWs 


(-2,-1) 


w 2 / 3 wlw 7 w 8 



(5.2) 



Higgsing b 

Model (b) corresponds to turning on vevs for H 23 , if 56 , H 78 and H^. The resulting 
brane tiling is shown in Figure |15 . 

The edge associated to H41 is contained in the original z\ path, which thus needs 
to be redefined. We can consider a new path z,\ given by Z\ = Z\Wi . The inte- 
grable system is summarized in fl5T3|) , where now (rii,n 2 ) corresponds to the (5™ 1 z^ 2 ) 
contributions. 
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1/4 0- 



7/8 



2/3 0—3- 



o- 



5/6 



1/4 0- 



Figure 15: Brane tiling for model (b) obtained by higgsing Y 4 ' . 



(ni,n 2 ) 


Loops 


(0,0) 


1 


(-1,0) 


W4/I + W 4 /lW 7 /s + W 2 /3W4/lW 7 / 8 + W 2 / 3 W4/lW 5 / 6 W 7 / 8 


(-2,0) 


W4/lW 7/8 + W 2 /3wl n W T /8 


(-2,1) 


W 4 /lW 7 / 8 


(-2,-1) 


W2/ 3 wl n W 7 /8 



(5.3) 



Higgsing c 

Finally, model (c) follows from turning on vevs for V15, V37, Vq 2 and V§4. The corre- 
sponding brane tiling is shown in Figure [Tj| 

This time both V15 and V37 overlap with the original z 2 path, which can be re- 
defined according to z 2 = z 2 wiw 3 . Denoting ''"■ " ■ ■ i1, << ' *" ! -" 2 
summarizes the resulting integrable system. 



rii,n 2 ) the (z? 1 ^ 2 ) contribution, (p~4| 



(ni,n 2 ) 


Loops 


(0,0) 


1 


(-1,0) 


W 8 /4 + W 3 / 7 W 8 /4 + W 6 / 2 W 3 / 7 W 8 /4 + 1 


(-2,0) 


W 8 /4 + W~fzW 8 /4 + W 3/7 W 8/ 4 + W~^W 3/7 W 8/ 4 + W~^ + W 3 /7™8/4 


(-3,0) 


W^W 8/ 4 + W~^W 3/7 W 8 /4 + ^3/7^8/4 + W l}h W V7 W l/A 


(-4,0) 


w-} b w 3/7 wl /A 


(-2,1) 


W 3 / 7 W 8 /4 



(5.4) 
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— Ch- 

3/7 

— •- 

— O- 
1/5 

— •- 



-4/8- 



— O- 

2/6 



— 6- 
-4/8- 



Figure 16: Brane tiling for model (c) obtained by higgsing Y 



4,0 



From an integrability point of view, this model is trivial, i.e. it only consists of 
Casimirs. From a quiver perspective, the reason for this is that, as one can deduce from 



Figure [T(| the associated quiver is fully non-chiral. This implies that all commutators 
vanish. 



6. Conclusions and Outlook 

We have investigated various applications of the correspondence between dimer mod- 
els and integrable system introduced by Goncharov and Kenyon in ||. We used it to 
explicitly construct relativistic generalizations of the periodic Toda chain associated 
to Y p, ° and Y p,p geometries. In these models, the calculation of commuting Hamilto- 
nians reduces to the combinatorics of non-intersecting paths on the brane tiling. We 
investigated the connection between the Kasteleyn matrix and the Lax operator, the 
non-relativistic limit of the integrable systems, additional relativistic versions of the 
periodic Toda chain based dimer models for general Y p,q geometries, and the identifi- 
cation of quiver impurities and spin chain impurities. Finally, we introduced a method 
for generating new integrable systems based on higgsing. We can envision, and are cur- 
rently pursuing, multiple directions in which the correspondence between dimer models 
and integrable systems can be exploited. We discuss some of them below. 

We have explained how higgsing is an efficient tool for generating new integrable 
systems. The characteristic polynomials for the dimer models associated to Z n x Z, m 
orbifolds of arbitrary geometries, which correspond to n x m arrays of copies of the 
original unit cell, can be determined using simple formulas ]l(J . These expressions have 
been used for calculating the multiplicity of perfect matchings associated to points in the 
toric diagrams of orbifolds p3| . It would be interesting to investigate whether there are 
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analogous expressions for the integrable systems associated to orbifolds starting from 
the integrable system for the unorbifolded geometry. New integrable systems could 
then be generated by higgsing, using the orbifold theories as starting points. 

It would be interesting to study the continuous limit of integrable systems arising 
from dimer models. By this, we not only mean the infinite length limit of integrable 
chains but also the (l+l)-dimensional field theory limit of fixed length systems. It 
is reasonable to conjecture that field theories such as A n Toda field theories can be 
constructed in this way. Furthermore, it is natural to expect that dimer models are 
useful for classifying (0+l)-dimensional, integrability-preserving defects and interfaces 
that can be added to such field theories. 



The work of [33| investigated the quantization of Riemann surfaces defined by 
the vanishing of the A-polynomials of three-manifolds M that are the complement 
of (thickened) knots or links. In the case of knots, the boundary of M is a 2-torus. 
The decomposition of M into glued tetrahedra with truncated vertices gives rise to 
a triangulation of the 2-torus, the developing map, whose dual is reminiscent of a 
dimer model. On the other hand, |5| and || discussed the relation between dimer 
models and quantum Teichmuller space. It is then natural to ask whether the existing 
similarities indicate the existence of a true connection that associates dimer models to 
three-manifolds that are knot complements. 

One clear direction for further research is to explore the connection between quan- 
tum integrable systems and gauge theories proposed by Nekrasov and Shatashvili @. 
By considering the 5d gauge theory on an ^-background with t\ — e and 62 = 0, the 
classical integrable systems we have investigated become quantized. We now review 
the interpretation of the quantization in terms of a B-brane on the spectral curve using 
the refined topological string pi. The mirror Calabi-Yau geometry takes the form 



uv + H(x,p) = 0. 

When the f2-background is turned on, the classical spectral curve, S defined by H (x, p) = 
0, is promoted to a quantum spectral curve 

H(x,p)^(x) = 

where \l/(x) is the wave-function for a B-brane. The quantum spectral curve is the 
Baxter equation for the relativistic Toda chain |2^|. The coefficients of H(x,p) are 
the energy eigenvalues of the Hamiltonians (and Casimirs) of the quantum integrable 
system. In the relativistic case, the differential operators e^ become shift operators, 
acting on \&(x) by 

e p ^(x) ->■ ^/(x + h). 
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Solutions to this differential equation were obtained from the refined topological string 
partition function in [J34fl . Thus the refined topological string should provide a way to 



solve the Baxter equation for quantum integrable systems. We plan to elucidate this 
connection in future work. 

A graphical representation for the finite non-periodic Toda lattice, similar to dimer 
models, was developed [^3| using the Poisson geometry of planar directed networks in 



an annulus fl36fl . The geometry of planar directed networks was developed in order 



to study the totally nonnegative Grassmannian [fj^] . Our work suggests an intriguing 
connection between the geometry of brane tilings and the totally nonnegative part of 
the double loop Grassmannian. 

In general, gauge theories arising from brane tilings have sequences of periodic 
Seiberg dualities known as duality cascades P7|-[f£T|. According to |5|], Seiberg duali- 



ties correspond to canonical transformations of the integrable system. In terms of the 
new variables the Hamiltonians will typically take a different functional form. However, 
since cascades are periodic, we are interested in special canonical transformations that 
preserve the functional form of the commuting Hamiltonians. Such canonical trans- 
formations are known as auto-Backlund-Darboux transformations. For example, the 
transformation ||3] 



_, _ dj + d-i ~ __ dj + d-i 

is an auto-Backlund-Darboux transformation of the relativistic Toda chain because it 
is canonical and the new Hamiltonian 



H 



^2 \ Ci + dr - 



takes the same functional form as the original Hamiltonian. The theory of auto- 
Backlund-Darboux transformations is closely related to the theory of separation of 
variables and discrete-time integrable systems [|3], f|5| . Thus we expect a fruitful inter- 



play between duality cascades and integrable systems. 
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